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We  report  first-principles  calculations  of  the  electronic  and  geometric  structure  of  the  (110)  cross-sectional 
surfaces  on  InAs/GaSb  superlattices,  and  compare  the  results  with  scanning  tunneling  microscopy  images  of 
filled  electronic  states.  In  both  the  predicted  and  measured  images  the  InAs  surfaces  appear  lower  than  GaSb, 
a  height  difference  we  show  is  caused  primarily  by  differences  in  the  electronic  structure  of  the  two  materials. 
In  contrast,  local  variations  in  the  apparent  height  of  surface  atoms  at  InSb-  or  GaAs-like  interfaces  arise 
primarily  from  geometric  distortions  associated  with  local  differences  in  bond  length.  Our  calculated  energies 
for  atomic  intermixing  indicate  that  anion  exchanges  are  exothermic  at  GaAs  interfaces  but  endothermic  at 
InSb  interfaces.  This  difference  may  explain  why  GaAs  interfaces  are  typically  more  disordered  than  InSb 
interfaces  in  these  heterostructures. 
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Cross-sectional  scanning  tunneling  microscopy  (XSTM) 
has  emerged  as  a  powerful  technique  to  characterize  III-V 
semiconductor  heterostructures. 1-11  Precise  characterization 
of  these  materials  is  made  possible  by  the  fact  that  a  zinc- 
blende  III-V  crystal  readily  cleaves  along  the  {110}  faces, 
producing  a  nearly  defect-free  surface  that  presents  a  cross- 
sectional  view  through  a  single  lattice  plane  of  structures 
grown  on  (001)  substrates.1  Tunneling  microscopy  is  particu¬ 
larly  useful  for  III-V  (110)  surfaces  because  of  the  simple 
surface  structure,  illustrated  in  Fig.  1(a).  The  III  atoms  relax 
towards  the  surface  and  V  atoms  away,  shifting  charge  be¬ 
tween  the  atoms  and  leaving  the  III  dangling  bond  essentially 
empty  and  the  V  surface  orbital  filled.  Because  the  STM 
surface  topography  in  constant-current  images  approximately 
corresponds  to  contours  of  constant  integrated  charge  den¬ 
sity,  only  the  III  dangling  bonds  are  seen  in  empty-state  im¬ 
ages  of  III-V  {110}  surfaces,  while  V  orbitals  are  seen  in 
filled-state  images.1212  Therefore,  XSTM  images  provide  an 
apparently  straightforward  chemical  identification  of  the  at¬ 
oms  observed. 

Since  the  first  report  of  atom-selective  STM  images  of 
GaAs(llO),12  and  the  observation  of  a  heterostructure  using 
XSTM,1  a  major  issue  has  been  delineating  between  elec¬ 
tronic  and  geometric  sources  of  height  contrast.  For  nomi¬ 
nally  homogeneous  materials  where  isolated  impurity  atoms 
are  observed,  such  as  dopants  or  substitutional  defects,  elec¬ 
tronic  origins  of  contrast  have  dominated  the  discussions.14 
For  heterostructures,  there  are  three  contrast  issues  to  be  con¬ 
sidered.  First,  the  different  III-V  materials  in  a  heterostruc¬ 
ture  usually  have  a  different  topographic  height  in  filled-state 
images.  Until  the  past  few  years, 6-10  discussion  of  this  dif¬ 
ference  focused  on  electronic  effects,  specifically  on  the 
band  gaps  and  band  alignments  (for  filled  states,  the  valence- 
band  maximum),  and  the  associated  number  of  bands  con¬ 
tributing  to  the  tunneling. 1-5  The  second  contrast  issue  is 
related  to  the  relative  appearance  of  point  defects  associated 
with  interdiffusion  between  the  materials.  For  example. 
Harper  et  al.  originally  described  As  defects  in  GaSb  as  ap¬ 
pearing  lower  in  height  because  of  the  position  of  the  As 


highest  occupied  molecular  orbital.15  Finally,  there  is  the  lo¬ 
cal  height  of  interfacial  bonds  to  consider.  For  systems  with¬ 
out  a  common  anion,  such  as  InAs/GaSb,  two  different  types 
of  interfacial  bonds  are  possible  (InSb  and  GaAs  bonds,  in 
this  case),  and  it  has  recently  been  proposed,  based  on  crys¬ 
tallographic  arguments,  that  the  local  XSTM  height  is  prima¬ 
rily  determined  by  local  bond  lengths.7'8 
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interface  interface  interface 


FIG.  1.  (a)  Relaxed  geometry  of  a  III-V  (110)  surface.  d0  de¬ 
notes  the  bond  length  in  the  bulk,  dt  the  out-of-plane  III-V  bond 
length  indicated,  and  the  height  difference  between  III  and  V 
surface  atoms,  (b)  Constant-current,  filled-state  XSTM  image  of  an 
InAs/GaAs  superlattice. 
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In  this  work,  we  use  first-principles  methods  to  describe 
the  electronic  and  geometric  structure  of  the  (110)  surfaces 
of  InAs/GaSb  superlattices.  We  find  that  the  apparent  surface 
height  difference  between  the  two  materials  is  primarily  an 
electronic  structure  effect,  but  the  local  height  differences 
observed  for  InSb  and  GaAs  interfacial  bonds  are  mostly 
geometric  in  nature.  In  addition,  the  calculations  reveal  that 
atomic  intermixing  lowers  the  energy  of  GaAs  interfaces, 
favoring  disorder. 

The  XSTM  measurements  were  performed  in  ultrahigh 
vacuum  using  InAs/GaSb  superlattice  samples  grown  by 
solid-source  molecular-beam  epitaxy,  as  described  in  detail 
previously.8  All  images  presented  here  are  of  (110)  surfaces 
recorded  with  filled  states  at  contant  current  (2. 2-2. 5  V, 
150-200  pA).  Figure  1(b)  shows  a  typical  filled-state  XSTM 
image  for  a  (110)  surface.  In  our  theoretical  work,  we  focus 
on  three  aspects  of  this  representative  image.  First,  the  large- 
scale  topography  shows  an  alternating  pattern  of  brighter  and 
darker  bands  (higher  and  lower  apparent  heights)  corre¬ 
sponding  to  GaSb  and  InAs  layers,  respectively.  Second,  a 
higher  Sb  row  is  evident  at  InSb  interfaces  in  this  figure, 
whereas  a  lower  As  row  is  evident  at  the  GaAs  interface. 
Third,  the  degree  of  atomic  disorder  at  the  two  interfaces  is 
qualitatively  different:  InSb  interfaces  typically  appear 
atomically  abrupt,  in  contrast  to  GaAs  interfaces  which  often 
exhibit  significant  disorder. 

Our  calculations  are  based  on  density-functional  theory, 
using  ultrasoft  pseudopotentials  as  implemented  in  the  VASP 
code.1617  Exchange-correlation  effects  were  treated  within 
the  local-density  approximation.  The  plane-wave  cutoff  was 
150  eV,  and  the  sampling  of  the  quasi-one-dimensional  Bril- 
louin  zone  was  equivalent  to  using  64  k-points  in  the  full 
zone  of  the  primitive  fee  cell. 

We  used  supercell  geometries  to  represent  the  InAs/GaSb 
heterostructures.  Because  of  the  importance  of  interfacial 
strain  in  this  material  system,  we  constructed  supercells  with 
starting  geometries  that  minimized,  as  much  as  possible,  any 
artificial  strain  at  the  interfaces.  Our  procedure  consisted  of 
the  following  three  steps. 

(1)  Bulk  calculations  were  performed  to  obtain  the  opti¬ 
mized  lattice  parameters  for  four  different  types  of  zinc- 
blende  BI-V  semiconductor  crystals:  InAs,  GaSb,  InSb,  and 
GaAs. 

(2)  For  each  of  these  four  different  homogeneous  materi¬ 
als,  we  constructed  slab  supercells  representing  the  unrecon¬ 
structed  (110)  surfaces.  This  was  done  by  periodically  repli¬ 
cating  the  unit  cells  obtained  in  the  previous  step  four  times 
along  the  (110)  direction,  resulting  in  slabs  containing  eight 
atomic  layers.  We  separated  adjacent  slabs  by  a  vacuum  re¬ 
gion  corresponding  to  five  atomic  layers,  which  we  con¬ 
firmed  was  sufficient  to  make  the  interaction  between  slabs 
negligible.  All  of  the  atoms  within  each  slab  were  then  re¬ 
laxed  within  the  constraint  of  the  fixed  in-plane  lattice  con¬ 
stant  determined  from  the  previous  step. 

(3)  Next  we  periodically  replicated  the  relaxed  slabs  six 
times  along  the  (001)  direction,  and  joined  two  such  ex¬ 
tended  slabs  together  to  form  various  III-V  heterostructures 
with  (001)  interfaces  and  exposed  (110)  surfaces.  For  each 
interfacial  bond  case  (InSb  and  GaAs),  the  two  slabs  were 


PHYSICAL  REVIEW  B  67,  121306(R)  (2003) 


TABLE  I.  Lattice  constants  a0  and  bond  lengths  dj  of  the  rel¬ 
evant  III-V  materials,  in  angstrom.  See  Fig.  1(a)  for  definitions  of 
the  different  bond  lengths. 


ao 

a0 

d0 

d\ 

d2 

InAs 

6.01 

6.06 

2.60 

2.62 

0.77 

GaSb 

6.04 

6.10 

2.62 

2.65 

0.76 

InSb 

6.43 

6.47 

2.78 

2.80 

0.86 

GaAs 

5.60 

5.65 

2.42 

2.44 

0.70 

“Lattice  constant  from  this  work. 
Experimental  lattice  constant. 


joined  at  a  distance  chosen  as  to  allow  every  interfacial  bond 
to  have  a  bond  length  corresponding  to  the  bulk  lattice  con¬ 
stant  computed  in  step  (1).  The  resulting  supercells  each  con¬ 
tained  a  total  of  96  atoms,  with  four  different  atomic  species, 
corresponding  to  a  (001)  superlattice  period  of  24  atomic 
layers.  Finally,  the  positions  of  all  atoms  were  completely 
relaxed  within  the  constraints  of  fixed  superlattice  period  and 
fixed  lattice  constant  along  [llO]. 

The  equilibrium  lattice  constants  calculated  in  step  (1)  are 
listed,  along  with  their  experimental  values,  in  the  first  two 
columns  of  Table  I.  The  agreement  is  very  good,  with  all 
errors  less  than  1%.  The  relaxed  (110)  surfaces  of  the  four 
homogeneous  materials,  as  obtained  from  step  (2),  show  the 
surface  buckling  obtained  in  many  previous  studies.13  The 
calculated  buckling,  illustrated  in  Fig.  1(a)  and  tabulated  in 
the  last  three  columns  of  the  table,  is  in  good  quantitative 
agreement  with  experiment;  in  particular,  the  height  differ¬ 
ence  between  III  and  V  atoms  for  InAs  is  0.77  A,  in  excellent 
agreement  with  the  value  of  0.78  A  determined  by  low- 
energy  electron  diffraction.18  After  relaxing  the  various  III-V 
heterostructures  described  in  step  (3),  we  simulated  XSTM 
images  using  the  method  of  Tersoff  and  Hamann.19  To  simu¬ 
late  filled-state  images,  we  integrated  the  local  density  of 
states  (LDOS)  from  1  eV  below  the  Fermi  level  up  to  the 
Fermi  level;  the  surface  of  constant  integrated  LDOS  then 
corresponds  to  the  ideal  STM  topography. 

Our  results  for  InAs/GaSb  heterostructures  with  InSb  in¬ 
terfacial  bonds  are  shown  in  Fig.  2.  The  geometry  of  the 
fully  relaxed  (110)  surface  is  displayed  in  Fig.  2(b).  As  in  the 
case  of  homogeneous  structures,  the  surface  atoms  buckle, 
causing  the  group-V  atoms  (As  and  Sb)  to  move  outward  and 
the  group-III  atoms  (In  and  Ga)  to  move  inward.  At  the  in¬ 
terface,  Sb  atoms  relax  still  further  outward  so  as  to  partially 
relieve  compressive  strain  in  the  InSb  bonds.  The  resulting 
simulated  XSTM  image  closely  resembles  the  measured  im¬ 
age,  as  indicated  by  the  inset  of  Fig.  2(a)  and  the  calculated 
XSTM  profile  across  the  row  maxima  shown  in  Fig.  2(c). 
Away  from  the  interface,  the  topographic  maxima  (from  the 
integrated  LDOS)  are  0.15  A  higher  on  the  GaSb  than  on  the 
InAs,  in  good  agreement  with  the  height  difference  of  about 
0.2  A  typically  observed  in  XSTM  images.  Interestingly,  the 
difference  in  height  between  the  actual  Sb  and  As  atoms 
associated  with  the  topography  is  much  smaller,  0.06  A, 
demonstrating  that  the  XSTM  height  difference  is  primarily 
caused  by  the  surface  electronic  structure.  Based  on  the  cal¬ 
culated  structure,  the  opposite  appears  true  for  the  local  to- 
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FIG.  2.  (a)  An  XSTM  image  of  an  InAs-GaSb  interface  with 
InSb  interfacial  bonds.  Inset:  Simulated  XSTM  image  of  this  struc¬ 
ture.  Both  gray  scales  span  about  1  A.  (b)  Side  view  of  the  fully 
relaxed  surface  geometry  for  this  interface,  (c)  Line  profile  for  the 
simulated  image  along  the  [001]  direction  (across  the  row  maxima). 
Circles  denote  the  relaxed  positions  of  the  surface  As  and  Sb  atoms, 
(d)  Line  profile  for  a  simulated  image  of  two  homogeneous  struc¬ 
tures  joined  together  at  their  ideal  lattice  positions  and  allowed  to 
relax  the  electronic  but  not  the  geometric  structure. 

pographic  height  difference  observed  at  the  InSb  interface. 
The  Sb  row  forming  InSb  bonds  is  geometrically  higher  by 
about  0.2  A  than  the  Sb  atoms  on  the  GaSb  surface,  essen¬ 
tially  the  same  height  difference  that  occurs  in  the  integrated 
LDOS,  demonstrating  that  the  observed  height  difference  in 
this  case  is  associated  with  the  geometric  structure. 

In  order  to  more  explicitly  delineate  the  relative  contribu¬ 
tions  of  geometric  and  electronic  structure  relaxation  on  the 
apparent  STM  topography,  we  calculated  the  electronically- 
relaxed  structure  of  an  “ideal”  InAs/GaSb  heterostructure, 
with  all  interfacial  atoms  frozen  at  ideal  positions.  In  this 
ideal  geometry  all  atoms  in  the  top  layer  have  exactly  the 
same  height,  and  thus  differences  in  heights  across  the  com¬ 
puted  topography  originate  from  purely  electronic  effects.  As 
shown  in  Fig.  2(d),  the  electronic  structure  alone  creates  a 
difference  in  height  of  0.13  A  between  the  InAs  and  GaSb 
surfaces,  close  to  the  0.15  A  difference  calculated  with  full 
relaxation,  further  supporting  our  conclusion  that  electronic 
structure  underlies  the  measured  height  difference  between 
the  two  materials. 

Our  analogous  results  for  InAs/GaSb  heterostructures 
with  GaAs  interfacial  bonds  are  summarized  in  Fig.  3,  with 
the  calculated  topography  again  in  qualitatively  good  agree¬ 
ment  with  that  observed  experimentally.  At  this  interface  the 
structural  relaxation  is  quite  different  from  the  InSb  case.  As 
shown  in  Table  I,  GaAs  bonds  have  the  shortest  surface  bond 
length  dx  and  thus  As  atoms  near  the  interface  relax  even 


FIG.  3.  (a)  An  XSTM  image  of  an  InAs/GaSb  interface  with 
GaAs  interfacial  bonds.  Inset:  Simulated  image  of  this  structure,  (b) 
Side  view  of  the  fully  relaxed  surface  geometry  for  this  interface, 
(c)  Line  profile  for  the  simulated  image  along  the  [001]  direction. 


further  inward  than  Sb  relaxes  outward  at  InSb  interfaces. 
The  resulting  geometric  height  of  the  As  row  at  the  GaAs 
interface  is  0.24  A  lower  than  the  As  atoms  far  from  the 
interface.  However,  as  we  saw  at  the  InSb  interface,  the  com¬ 
puted  topographic  height  difference  is  almost  identical,  0.23 
A,  indicating  that  the  local  depression  of  the  GaAs  interface 
is  almost  completely  geometric  in  origin. 

Finally,  we  address  why  interfacial  roughness  appears  to 
depend  on  the  interfacial  bond  type,  with  GaAs  interfaces 
generally  observed  to  be  more  disordered.  Experimentally, 
most  defects  occur  close  to  the  interfaces,  suggesting  that 
they  arise  from  simple  Ga-In  or  As-Sb  exchanges  across  the 
interface,  rather  than  from  bulk  defects  such  as  vacancies  or 
cation-anion  antisites.  An  example  of  an  apparent  Sb  atom 
observed  in  an  As  site  at  a  GaAs  interface  is  highlighted  in 
Fig.  4(a).  To  confirm  the  structural  assignment  of  such  fea¬ 
tures,  we  theoretically  modeled  such  a  defect  by  replacing 
one  of  the  surface  As  atoms  at  a  GaAs  interface  with  an  Sb 


FIG.  4.  (a)  An  XSTM  image  of  InAs/GaSb  with  a  GaAs  inter¬ 
face.  An  apparent  Sb  atom  in  an  As  site  at  the  interface  is  circled. 
Inset:  Simulated  XSTM  image  for  such  a  defect,  (b)  Top  view  of  the 
Sb-in-As-site  defect  model  structure,  with  the  top-layer  atoms 
shown  larger,  (c)  Model  used  to  calculate  the  energetics  of  an  As-Sb 
exchange  across  an  InSb  interface,  (d)  Model  for  an  As-Sb  ex¬ 
change  across  a  GaAs  interface. 
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atom,  as  illustrated  in  Fig.  4(b).  The  simulated  XSTM  image 
for  this  structure,  shown  in  the  inset  of  Fig.  4(a),  closely 
resembles  the  experimental  result. 

To  investigate  the  energetics  of  anion  interfacial  defects, 
we  consider  the  simplest  defects  that  both  preserve  the  global 
stoichiometry  and  satisfy  local  chemical  bonding  require¬ 
ments,  As-for-Sb  exchanges.  Such  exchanges  represent  a 
simple  mechanism  for  interfacial  disorder  at  a  nominally 
abrupt  interface.  We  studied  the  energetics  of  exchanging 
adjacent  As  and  Sb  atoms  both  at  an  InSb  interface  and  at  a 
GaAs  interface,  as  shown  in  Figs.  4(c)  and  4(d),  respectively. 
The  structures  were  fully  relaxed  before  and  after  the  ex¬ 
change  and  the  change  in  the  total  energy  then  computed.  At 
the  InSb  interface,  the  exchange  raised  the  total  energy  by  7 
meV.  Surprisingly,  at  the  GaAs  interface,  the  same  process 
actually  lowered  the  total  energy  by  22  meV;  i.e.,  the  forma¬ 
tion  of  such  defect  pairs  is  exothermic.  This  result  implies 
that  abrupt  GaAs  interfaces  are  thermodynamically  unstable. 
Therefore,  although  kinetic  barriers  may  suppress  anion  ex- 


PHYSICAL  REVIEW  B  67,  121306(R)  (2003) 

changes,  one  should  generally  expect  GaAs  interfaces  to  be 
more  disorderd  than  InSb  interfaces  (as  widely  observed). 

In  summary,  we  have  used  first-principles  electronic- 
structure  methods  to  clarify  the  interpretation  of  XSTM  im¬ 
ages  of  (110)  surfaces  on  cleaved  InAs/GaSb  heterostruc¬ 
tures,  focusing  on  the  differences  between  interfaces  with 
InSb  versus  GaAs  bonds.  We  find  that  the  apparent  height 
differences  between  the  InAs  and  GaSb  surfaces  are  largely 
associated  with  the  electronic  structure,  whereas  the  local 
height  differences  at  the  InSb  and  GaAs  interfaces  are  caused 
by  geometric  relaxation  from  the  partial  relief  of  local  bond 
strain.  Finally,  our  calculations  of  the  energies  associated 
with  interfacial  exchange  of  anions  reveal  that  As-for-Sb  ex¬ 
changes  are  exothermic  at  GaAs  interfaces,  but  endothermic 
at  InSb  interfaces.  This  result  may  explain  why  GaAs  inter¬ 
faces  are  generally  more  disordered  than  InSb  interfaces  in 
InAs/GaSb  heterostructures. 
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